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SYSTEM AND METHOD FOR DETERMINING A FLOW PROFILE IN A DEVIATED 

INJECTION WELL 

CROSS-REFERENCE TO RELATED APPLICATIONS 

This application claims priority from US Provisional Application 60/510,596, filed 
October 10, 2003, which is incorporated herein by reference. 

BACKGROUND OF THE INVENTION 

Field of the Invention 

[0001] The present invention relates to a system and methodology for determining 

a flow profile in a well, and particularly to determining a flow profile in a deviated 
injection well. 

Description of Related Art 

[0002] In a variety of wells, various parameters are measured to determine specific 

well characteristics. For example, temperature logging has been used for profiling the 
injection rate in vertical wells. Existing methods of analyzing injection profiles are 
designed for vertical wells where the injection interval is usually small and the time to 
flush the wellbore volume is negligible. Also, the displacement process of the reservoir 
fluid can be represented by a radial flow model. 

[0003] However, if the wellbore is deviated, such methods do not enable profiling 

of the injection rate. Thus, deviated wellbores, such as horizontal wellbores, present 
greater problems in evaluating and predicting flow profiles for injection wells. 
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BRIEF SUMMARY OF THE INVENTION 

[0004] In general, the present invention provides a system and methodology for 

using a well model in determining characteristics of an injection well. For example, the 
system and methodology enables the use of temperature profiles in a deviated injection 
well to determine a flow profile in such well. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0005] Certain embodiments of the invention will hereafter be described with 

reference to the accompanying drawings, wherein like reference numerals denote like 
elements, and: 

[0006] Figure 1 is an elevation view of a completion and sensing system deployed 

in a deviated wellbore, according to an embodiment of the present invention; 

[0007] Figure 2 is an elevation view of the system illustrated in Figure 1 

illustrating a fluid being injected into the well, according to an embodiment of the present 
invention; 

[0008] Figure 3 is a flowchart generally representing an embodiment of the 

methodology used in determining a flow profile in a well, according to an embodiment of 
the present invention; 

[0009] Figure 4 is a diagrammatic representation of a processor-based control 

system that can be used to carry out all or part of the methodology for determining flow 
profile in a given well, according to an embodiment of the present invention; 
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[0010] Figure 5 is a flowchart generally representing a methodology for 

determining flow profiles based on temperature profiles during fluid injection into a 
deviated well, according to an embodiment of the present invention; 

[0011] Figure 6 is a graphical representation plotting temperature against distance 

along a wellbore during an early period of injection; 

[0012] Figure 7 is a graphical representation similar to that of Figure 6 but with a 

different injection geometry; 

[0013] Figure 8 is a graphical representation similar to that of Figure 6 but at a 

later injection time; 

[0014] Figure 9 is a graphical representation similar to that of Figure 7 but at a 

later injection time; 

[0015] Figure 10 is a flowchart generally representing a methodology for 

determining flow profiles based on temperature profiles in a deviated well during a shut-in 
period, according to an embodiment of the present invention; 

[0016] Figure 1 1 is a graphical representation plotting temperature against distance 

along a wellbore during a shut-in period; 

[0017] Figure 12 is a graphical representation similar to that of Figure 1 1 but with 

a different injection geometry; 

[0018] Figure 13 is a schematic representation of a processor system that receives 

data related to temperature profiles and other well parameters to derive flow profiles; 

[0019] Figure 1 4 is a schematic representation of a deviated well divided into a 

multi-segment grid system for modeling; 
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[0020] Figure 15 is a graphical representation of dimensionless temperature 

plotted against dimensionless time; and 

[0021] Figure 16 is a flowchart generally representing a methodology for 

determining flow rates for a plurality of intervals along a deviated well. 

DETAILED DESCRIPTION OF THE INVENTION 

[0022] In the following description, numerous details are set forth to provide an 

understanding of the present invention. However, it will be understood by those of 
ordinary skill in the art that the present invention may be practiced without these details 
and that numerous variations or modifications from the described embodiments may be 
possible. 

[0023] The present invention generally relates to a system and method for 

determining flow profiles in a deviated well. Temperature measurements are taken along 
a wellbore, and those measurements are vised in determining flow profiles along a deviated 
injection well, such as a generally horizontal injection well. In some applications, a flow 
profile is derived based on data obtained during injection of a fluid into the deviated well. 
In other applications, a flow profile is derived based on data obtained during a shut-in 
period following injection or during periods of resumed injection. 

[0024] A temperature sensing system, such as a distributed temperature sensor, is 

deployed with an operational completion and enables temperature measurements to be 
taken during fluid injection periods or during shut-in periods. Based on the collected 
temperature data, flow profiles of the injected fluid along the deviated well can be derived. 

[0025] In general, when a cool fluid, such as a liquid, e.g. water or oil, or a gas, is 

injected into a hot reservoir, a variety of thermal changes occur. For example, during 
injection, cool fluid moves through the wellbore and into the reservoir while heat flows 
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from the reservoir toward the wellbore. A similar effect occurs along the axis of the 
wellbore as fluid flows from the heel of the wellbore toward the toe, and heat flows from 
the toe of the wellbore toward the heel. The thermal characteristics of the heat flow can be 
modeled in a manner that enables determination of the flow profile of the fluid flow into 
the reservoir. Other factors, such as thermal conductivity of the surrounding formation, 
may also be utilized in modeling the flow profile, as discussed below. 

[0026] Furthermore, if the injection of cool fluid is stopped, thereby creating a 

shut-in period, other unique thermal changes occur that enable the determination of 
injection profiles. For example, when the injection stops, the wellbore begins heating up, 
but not necessarily uniformly. The temperature recovery provides an indication of where 
the cool fluid was moving into the reservoir during injection. Reservoir intervals that 
were receiving substantial amounts of injection fluid, for example, are slower to rise in 
temperature during the shut-in period. These thermal changes are applied to models that 
enable the derivation of flow profiles, as discussed below. 

[0027] Referring generally to Figure 1 , a system 20 is illustrated in accordance 

with an embodiment of the present invention. System 20 comprises a completion 22 
deployed in a well 24. In this example, well 24 is a deviated well having a generally 
vertical section 26 and a deviated section 28, such as the generally horizontal section 
illustrated in Figure 1 . Well 24 is defined by a wellbore 30 drilled in a formation 32 
having, for example, one or more fluids, such as oil and water. A tubing 33 extends 
downwardly into wellbore 30 from a wellhead 34 disposed, for example, along a seabed 
floor or a surface of the earth 36. In the illustrated example, tubing 33 extends to a casing 
shoe 38 that may be located at the lower end of vertical section 26 above a heel 40 of 
deviated section 28. Completion 22 is disposed in deviated section 28 and may extend 
from casing shoe 38 through heel 40 toward a toe 42 of well 24. In many applications, 
wellbore 30 is lined with a casing 44 that may be perforated to enable fluid flow 
therethrough. 
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[0028] As further illustrated, system 20 comprises a temperature sensing system 

46. For example, temperature sensing system 46 may comprise a distributed temperature 
sensor (DTS) 48 that is able to continually sense temperature along deviated section 28 of 
wellbore 30 at multiple locations. Distributed temperature sensor 48 may be coupled to a 
controller 50 able to receive and process the temperature data obtained along wellbore 30. 
As discussed in greater detail below, controller 50 also enables use of the temperature data 
in conjunction with a model of the well to derive injection flow profiles of fluid flowing 
from completion 22 into formation 32 along deviated section 28 of the well 24. 

[0029] During injection, a fluid, such as water, is pumped down through tubing 33 

and into completion 22 along the illustrated horizontal section of the well. The fluid is 
forced outwardly along deviated section 28 such that fluid flows from completion 22 into 
formation 32, as indicated by arrows 52 of Figure 2. The injectivity profile or flow profile 
of fluid moving from the wellbore 30 into formation 32 often is not uniform. For 
example, the flow may be substantially greater in proximity to heel 40 as opposed to toe 
42. Also, the formation material may vary along the length of deviated section 28, which 
can extend a substantial distance, e.g. up to several kilometers. 

[0030] Referring generally to Figure 3, an example of the methodology of the 

present invention is illustrated in flow chart form. Determining flow profiles along a 
given injection well comprises deploying a sensor system in the well with an operable 
injection completion, as illustrated by block 54. Then, an injection fluid is injected into 
formation 32 via completion 22, as illustrated by block 56. The sensor system may 
comprise a distributed temperature sensor designed to sense temperature along deviated 
wellbore section 28, as illustrated by block 58. As discussed more fully below, the 
sensing of well parameters can be done during injection and/or subsequent to injection 
during a shut-in period. A well model may then be applied to determine the flow profile 
of injected fluid along deviated section 28 of well 24, as illustrated by block 60. 

[0031] Some or all of the methodology outlined with reference to Figures 1 -3 may 

be carried out by controller 50 which comprises an automated system 62, such as the 
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processing system diagramatically illustrated in Figure 4. Automated system 62 may be a 
computer-based system having a central processing unit (CPU) 64. CPU 64 may be 
operatively coupled to temperature sensing system 46, a memory 66, an input device 68, 
and an output device 70. Input device 68 may comprise a variety of devices, such as a 
keyboard, mouse, voice-recognition unit, touchscreen, other input devices, or 
combinations of such devices. Output device 70 may comprise a visual and/or audio 
output device, such as a monitor having a graphical user interface. Additionally, the 
processing may be done on a single device or multiple devices at the well location, away 
from the well location, or with some devices located at the well and other devices located 
remotely. 

[0032] In automatically determining a flow profile or profiles of fluid injected into 

formation 32, a model utilizing temperature changes along deviated section 28 as an 
indicator of flow profiles may be stored by automated system 62 in, for example, memory 
66. As illustrated best in Figure 5, the general approach involves obtaining an initial 
temperature profile along at least deviated well section 28, as indicated by block 72. 
Additionally, the temperature of the injected fluid, e.g. water, is measured, as indicated by 
block 74. This value also may be stored on automated system 62 for use in modeling the 
injectivity profile. In many applications of the modeling technique, a greater contrast 
between the temperature of the injected fluid and the temperature of the reservoir can 
improve the usefulness of the model. By way of example, the injected fluid may be at a 
temperature of 60-70 degrees Fahrenheit, and a reservoir may be at a temperature of 200- 
240 degrees Fahrenheit. However, the temperature of the injected fluid and the reservoir 
can vary substantially from one application and environment to another. The fluid is thai 
injected, as illustrated by block 76. Subsequently, temperature profiles are obtained along 
the deviated section 28, and this data is provided to automated system 62 via the 
temperature sensing system 46, as illustrated by block 78. The temperature profiles can be 
taken during injection or during a shut-in period subsequent to injection depending on the 
particular model applied. In either event, the model is used to derive a flow profile of the 
injected fluid, as illustrated by block 80. For example, the data collected may be 
processed according to the model/algorithm stored on automated system 62 to 
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automatically present a well operator with detailed information on the injection flow 
profile via, for example, output device 70. 

[0033] In determining the flow profiles, a given a well model may utilize thermal 

behavior characteristics that occur during injection. An example of a well in which 
injectivity decays along a horizontal well axis is illustrated in Figure 6 which provides a 
graph of temperature plotted against distance along the horizontal section of the wellbore. 
The graph illustrates temperature changes along the horizontal wellbore section during the 
first few hours of injection. As illustrated by a graph line 82 at 0.002 day, graph line 84 at 
0.01 day, graph line 86 at 0.03 day, graph line 88 at 0.05 day, graph line 90 at 0.07 day, 
graph line 92 at 0.1 day, and graph line 94 at 0.2 day, a gradual propagation of the 
temperature front in the wellbore is observed due to a much greater injectivity towards 
heel 40 than toe 42. A cooling or reduction of the temperature at heel 40 is also 
graphically illustrated. 

[0034] The thermal characteristics of a well in which injectivity is skewed towards 

toe 42 is illustrated in Figure 7. In this example a much faster propagation of the wellbore 
temperature front is observed, as indicated by the location of graph lines 82, 84, 86, 88, 
90, 92, and 94. 

[0035] Referring generally to Figures 8 and 9, temperature profiles for the 

examples provided in Figures 6 and 7, respectively, are again graphically illustrated, but at 
a substantially later period of the injection. In each of these figures, a graph line 96 
represents a temperature profile after one day of injection, and a graph line 98 represents a 
temperature profile after two days of injection. The thermal characteristics, such as those 
illustrated in Figures 6, 7, 8 and 9, demonstrate thermal changes that occur during 
injection of fluid into formation 32. Those thermal changes can be used by an appropriate 
model to determine flow profiles. In other words, the injection of a cooler fluid into the 
reservoir at different rates along deviated section 28 creates thermal changes over the 
period of injection. By accurately measuring the thermal changes via, for example, 
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distributed temperature sensor 48, the actual injection flow profiles can be derived by the 
appropriate model . 

[0036] At other phases of the process, useful thermal data also may be obtained. 

For example, following an injection phase, a shut-in phase can lead to interesting thermal 
events which can be modeled to provide an injection flow profile. When the injection of 
fluid starts, the wellbore begins heating, but not necessarily uniformly. The taking of 
temperature profiles during this temperature recovery period provides an indication of 
where the cooler injection fluid is moving into the reservoir dialing injection. For 
example, reservoir intervals receiving a greater flow of the cooler fluid are slower to 
regain heat during the shut-in period. Thus, the temperature profiles taken during a shut- 
in period can be used to determine injection flow profiles. 

[0037] As illustrated best in Figure 10, the general methodology for utilizing shut- 

in data involves obtaining an initial temperature profile along at least deviated well section 
28, as indicated by block 100. The temperature of the injected fluid, e.g. water, also is 
measured, as indicated by block 102. The injection fluid is then injected, as illustrated by 
block 104. Following an injection period, e.g. two days of injection, the injection is shut 
down for a shut-in period, as illustrated by block 1 06. Subsequently, temperature profiles 
are obtained along the deviated section 28, as illustrated by block 108. This data, along 
with other collected data, may be provided to automated system 62 via the temperature 
sensing system 46. By applying the appropriate model, the shut-in temperature data can 
be utilized in deriving an injection flow profile along the deviated section 28 of well 24. 
In some applications, the injection is resumed, and that resumed injection may be followed 
by a subsequent shut-in period, as indicated by block 110. The repeating of injection and 
shut-in periods can be used to obtain additional data, to verify results, and/or to 
continually monitor the injection flow profile. It should also be noted that appropriate 
models can also be designed to utilize the thermal characteristics of a well when injection 
is resumed after a shut-in period. As the fluid is reinjected, the well warms up and a step 
rise in temperature is indicated. This step or slug moves as a front along the deviated 
section of the wellbore and provides an indication of the flow profile. If, for example, the 
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front moves slowly, this generally indicates greater flow towards the heel of the deviated 
section. If, on the other hand, the front moves more rapidly, this can indicate greater flow 
toward the toe of the wellbore. 

[0038] In determining the flow profiles based on data obtained during the shut-in 

period, the well model utilizes thermal characteristics that occur during shut-in. In Figures 
1 1 and 12, graphs of temperature plotted against distance along the horizontal section of 
the wellbore are provided for the scenarios described above with reference to Figures 6 
and 7. However, the data graphically illustrated in Figures 1 1 and 12 represents 
temperature profiles taken during a shut-in period following the injection period illustrated 
graphically for a first scenario in Figures 6 and 8 and for a second scenario in Figures 7 
and 9, respectively. The graphs of Figures 1 1 and 12 illustrate temperature changes along 
the horizontal wellbore section at various time points of the shut-in. Specifically, the 
temperature changes are indicated by a graph line 112 providing a temperature profile at 
the start of the shut-in period, a graph line 114 providing a temperature profile at 0.5 days 
into the shut-in, and a graph line 116 providing a temperature profile at 1 day into the 
shut-in period. From the temperature profile data, it becomes apparent that the 
temperature rebounds quickly at intervals of the deviated wellbore having a lower rate of 
injection. On the contrary, intervals with greater injectivity, i.e. a greater rate of flow into 
the reservoir, rebound more slowly. The differences in thermal characteristics of the 
temperature recovery along the length of the deviated well enable determination of the 
injection flow profile. 

[0039] In some applications, the accuracy of the flow profiles can be improved by 

accounting for additional well related parameters. As illustrated in Figure 13, the use of a 
subject model 118 can include additional inputs other then the primary input of 
temperature profiles 120. In this example, the model is utilized or processed on automated 
processor system 62, and a variety of data is fed into the model and processor system 62 
via, for example, sensors or manual input via input device 68. For example, temperature 
profile data 120 may be provided by distributed temperature sensor 48. Other well related 
parameters, such as recent history 122, permeability of the reservoir 124, injection rate 
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1 26, injection period 1 28, and/or thermal conductivity 130, can be utilized by model 1 1 8 
on processor system 62 to provide reliable injection flow profiles to a well operator. 

[0040] A specific model/algorithm for determining flow profiles based on thermal 

data obtained during injection of fluid may take a variety of physical phenomena into 
account. For example, the injection of a cool fluid into a relatively hot reservoir creates 
both a flow of fluid and a flow of heat. Cool or cold fluid moves through the wellbore and 
into the reservoir as heat flows from the reservoir toward the wellbore. A similar effect 
occurs along the wellbore axis in that fluid flows from the heel to the toe, and heat flows 
from the toe to the heel. 

[0041] A practical model,to predict the temperature distribution along the wellbore 

when the injection flux is specified, or to estimate the injection flux distribution with 
measured temperature profile can be described as follows. Initial assumptions are set 
forth in the following points: 

(1) The pressure and temperature gradient in the direction parallel to the axis of the 
horizontal well (x) are much smaller than those in the directions (y and z) 
perpendicular to the axis. So, the mass and heat transference along x-direction can 
be neglected. 

(2) The viscosity and density of oil and water are constants. So, the pressure and water 
saturation distribution are independent of temperature. 

(3) Isotropic formation. 

(4) Incompressible injected fluid. 

(5) Open hole completion. 

[0042] The wellbore model in this embodiment further utilizes a wellbore flow 

rate distribution equation and a temperature distribution equation described below. 
(1) Wellbore flow rate distribution equation 

Let Q wi (x 9 t) denote the injected fluid volume flow rate along the 

wellbore (m 3 I hour), and q wi (x,t) denote the fluid volume rate 
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injected into the formation for unit length of the wellbore 
(m 2 / hour). The mass conservation of injected fluid in the 
wellbore gives: 

Q wi 0 - Q wi (x + dx,t)- g w . (x, t)dx = 0 



8x 



■ = -q«(.x t t) 



Integrating equation (1.1) from the heel to toe yields: 



Q bv (t) = Q wi (0,t) = jq wi (x,t)dx 

0 

Q wl {x,t) = Q iaj {t)-\q wl ^,t)^ 



where Q M (t) is the total injection rate at the heel of a horizontal well. 



(1.1) 



(1.2) 




Q v (x,t) 



Wellbore flow rate distribution 



(2) Temperature distribution equation: 
The total heat stored in the tiny element dx is: 

Q well =c i A w T w (x,t)dx (1.3) 

Let q Tw (x,t) = 17U-A — I denote the heat flow rate from the formation into the unit 



\ dr 

length wellbore and T w (x,t) denote the temperature profile along the wellbore. The 
energy conservation equation for element dx is: 



dt 



*7 win ^ y 



(1.4) 
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where q^ in and denote the heat rate flowing into and out of the tiny wellbore element 

dx. q^ in is composed of two terms: the heat carried by the fluid flowing into the element 

through the wellbore cross-sectional area at x, c.Q^(xj)T w (xj), and the heat flowing 

from the formation to the wellbore element through the wellbore surface due to heat 
conduction, q Tw (x,t)dx . q wout is also composed of two terms: the heat carried by the fluid 

flowing out of the element through the wellbore cross-section area at jc + dbc, 
CiQ^ix + dXttyTwix + dxj), and the heat carried by the fluid flowing out of the element 
through the wellbore surface, c^.foOT^foOdx . 

=c i Q wi {x,t)T w {xj)-\-q Tw {xj)dx 



-(f).. 



?7v(*.0 = 2>n7r, 

= CiQ* (* + d*. i)T w O + dx, t) + c iqwi (x, t)T w (x, t)dx . 
Substituting q Mn , q ml and equation (1 .3) into (1 .4) gives: 

SiQ^TJ (BT\ . dT w 



dx y 8r J r=rw dt 

Substituting equation (1.1) into the above equation yields: 



where c t -- heat capacity of injected fluid (J/(m\°iO) ; 
A w — flowing area in the wellbore (m 2 ) ; and 
77 — thermal conductivity of the formation (J/(m.° K.h)) . 



[0043] Near Wellbore Heat Transportation - Water Injection 

The near wellbore flow regime can be regarded as steady-state radial flow. Consider a 
tiny radial element between r and r + dr . The total heat stored in this element for water 
injection is composed of three terms: the heat stored in the water phase 
Q w = 27rrdr(ps w (x,r,t)cjr(x,rj) , the heat stored in the oil phase 

Q o =27irdr0\l-s w (x,r,t))c o T(x,r,t) , and the heat stored in the rock 

Q r = 2nrdr(l-0)c r T(x,r 9 t) , i.e., 
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Q = Q +Q +Q 

= 27irdrT(x, r, t)[s w (x, r, 0*(c w - c Q ) + c> + (1 - ^)c r ] 

where T(x,r,t) is the temperature distribution in the reservoir, ^ is porosity, c Q is the 

heat capacity of oil (J/(m 3 .° K)) , and s w (x,r,t) is water saturation. 

3T 

Let q T (x,r,t) = 2x7} r denote the heat inward radial flow rate, q w (x^r,t) denote the 

dr 

water volume outward radial flow rate, and ? 0 (x,r,r) denote the oil volume outward 
radial flow rate. The energy conservation equation for this radial element is: 



dQres 



dt 

where q Hn and q rout denote the heat rate flowing into and out of the tiny radial element 
dr . q Hn is composed of two terms: the heat carried by the oil and water flowing into the 
element through the inner surface at r, [c w q w ( x , r J) + c 0 q 0 (x,r,tjft{x,r,i) , and the heat 
flowing into the element through the outer surface at r + dr due to heat conduction, 
q T (x, r + dr, t) . q rout is also composed of two terms: the heat carried by the oil and water 

flowing out of the radial element through the outer surface at r + dr, 
[c w # H ,(* 5 r + dr,0 + c o £ a (x,r + dr,t)]T(x,r + dr,t), and the heat flowing out of the tiny 
element through the inner surface at r due to heat conduction, q T (x, r, t) . 

<lrin = kvtfw (*> r * 0 + c o<lo (*> r > OFfo r,t) + q T (x,r + dr, 0 
9n*r = i c ^ w (x* r + dr > 0 + c * (*■ r + dr > 0F(*. r + dr, t) + ? r (r, f) 

g r (x, r,0 = 2^r— . 

or 

Substituting the above equations and equation (1.6) into equation (1.7) yields: 

(1.8) 



k (x, r, 0*(c w - c c ) + cj + (1 - <*)c r ]^ + *(c w - c Q )T(x, r, /) ^ 
1 



^[r? u%n -c„ 8( ^> 



2^r ' dr\^ dr ) w dr ° dr 
Let r w/ (x, 0 denote the water front at time / , and choose r < r w/ (x, t) . Then: 

■s w (x,r, t) = 1 - s or , ? w (x, r, f ) = ^ w . (x, 0 , q o (x, r, t) = 0 
Near the wellbore, the heat flux also can be considered as a constant, i.e., 
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dr\ dr J dr 



Equation (1.8) becomes: 



[a-,^c,-o + ^ +( i-^]f = _£^o^ 

At the wellbore ( r = r w ), and: 
[d - s or Mc w -c 0 ) + cj + (l -<j>)c r f^ = - ^^ (f) • (1-10) 
Combining equation (1.10) with equation (1 .5) and choosing c. = c w yields: 
{ \ [(1 - 'a, MP W - c 0 ) + cj + (1 - flc, ] + c^J^f = ~^g w , ^ • OH) 



Let: 



Z4, i+a, w g 2 4;r77 



Equation (1.11) becomes: 

—?L + La w (x,t)—^ = 0 (1.12) 
3f fix 



[0044] Near Wellbore Heat Transportation - Oil Injection 

For oil injection, the water phase flow rate is q w {x,r,t) = 0 and water saturation is 

s w (x,r,t) = s wi . The total heat stored in the radial element is composed of three terms: the 
heat stored in the water phase Q w = 27irdr<f)s wi c w T(x,r,t) , the heat stored in the oil phase 
Q o = 27zrdr0[l-s wi )c o T(x,r,t) , and the heat stored in the rock 
Q r = 2nrdr(\-<f>)c r T{x,r,t) , i.e., 
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Q res =Q W +Q 0 +Q r 

= 27trdrT{x, r, t)[s wi (/>(c w - c Q ) + cj + (1 - 0)c r ] 



Similarly: 



<lrin = c o<lo ( x > r > 0T(x> r,t) + q T (x 9 r + dr, 0 
Rrout = c o q o (x, r + dr, t)T(x, r + dr,t) + q r (r, t) . 
The energy conservation equation for this radial element is: 
Substituting the above equations into equation (1 .7) yields: 



dt 27TT 

In the steady state flow regime near the wellbore: 



2nrj 



dr 



r 

\ dr J 



— c 



d(q 0 T) 
0 dr 



q o (x,r,t)=q wi (x,t), 



dr 



' drj 



= 0. 



And, equation (1.13) becomes: 



r j., ^ j. „ A idT c c q wi (x,t)dT 

M(c w -c 0 ) + cj + (l- <f>)c r J— = o ^: ,v — 

dt 2nr or 



At the wellbore ( r = r w ): 



Combining equation (1.15) with equation (1.5) and choosing c, = c a yields: 

f {i7<pL [s ^ Cw _ Ce ) + cJ + (1 _ , K ]+ - -c o& a7: 



Let 



C w C w A w 



<*„(*,*) = ■ 



LA.. 



\ + a ]o a 2 



Equation (1.16) becomes: 



—£- + La(x,t) — ^ = 0 
dt ° dx 



(1.13) 



(1.14) 



M(c w -cJ + C > + (l-^J^ = -^^[f) ; (1.15) 



(1.16) 



(1.17) 
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[0045] Boundary Condition and Initial Condition 

The boundary condition is specified with the injection temperature at the heel ( x = 0): 

^(0,0 = 7^(0 0-18) 
The temperature at the heel T^if) can be determined with a wellbore heat 

transmission model, such as the H. J. Ramey model. 
The initial condition is: 

T w (x,0) = T R (1.19) 

where T R is the reservoir temperature. 

[0046] Wellbore Temperature Prediction - Forward Problem 

When the reservoir properties (porosity (f> and permeability k\ fluid properties 

(density p o and p w , viscosity /j 0 and ju w , relative permeability k ro and ), wellbore 

geometry (wellbore diameter r w , fluid flow area A w , roughness s 9 and the length of 

perforated section L) are specified, then the injection rate distribution q wi (x,t) can be 

determined with an analytical model, such as the model established by TUPREP, or a 
numerical model, such as the ECLIPSE100. And thus, with the properly defined boundary 
condition (1.18) and initial equation (1.19), the wellbore temperature profile T w (jx 9 t) can 

be predicted by solving equation (1.12) for water injection or equation (1.17) for oil 
injection. 

Let q 0 = ^ inj ^ °^ denote the average injection flux at time t =t 0 , and 
L 

a a 

o w 



A w \l + a lw a 2 



C w#0 J 

Equations (1 .12) and (1.17) can be rewritten as: 



V 0 ( 4^ 



C w < lo ) 



^ + a Dw (x,t)La w ^ = 0 (1 .20) 

dt ox 
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^ + a Do (xj)La o ^ = 0 . (1.21) 
dt ox 

m 2 1 

The unit of q wi (x,t) and q 0 is , so the unit of a w and a is . 

hour hour 

x 

Let £• = —-, t D = <2 w ? for water injection or / D = aj for oil injection. They denote the 
L 

dimensionless variables. The dimensionless forms of equations (1.20) and (1.21 ) are: 

-f^ + a^t^-^^O (1.22) 
dt D dg 

—^ + a Do (g,t D )—^ = 0 . (1.23) 
dt D dg 

Or, both equations can be rewritten as: 

^L + a Di (g,t D )^ = 0; (i = o,w) (1.24) 
dt D dg 

Let g c {t D ) denote the characteristic curve along which the temperature is unchanging, 

i.e., 

dT w =-^ + ^^1 = 0 . (1.25) 
dt D dg c dt D 



Comparing equation (1 .24) with (1 .25 ) yields 



= ci Di {g c ,t D ) (1.26) 



Equation (1.26) is the characteristic equation with respect to the partial differential 
equation (1 .24). Equation (1 .26) defines a group of curves, characteristic curves. It can be 
proved that all characteristic curves do not intersect with each other. If one characteristic 
curve crosses the positive g coordinate, then the temperature on this curve is specified by 

the initial condition, i.e., equal to the reservoir temperature T R . Otherwise, the curve will 

cross the positive t D coordinate, and the temperature on this curve is specified by the 

boundary condition T w (t Dp ) , where t Dp is the intersection of the characteristic curve with 

the time coordinate. 
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[0047] The modeling technique described above enables the determination of 

injection flow profiles based in large part on temperature profiles obtained during 
injection of the fluid. However, the shut-in period also can be modeled such that injection 
flow profiles can be determined based on thermal information obtained during the shut-in 
period. Of course, the data obtained and modeled during the injection period and the shut- 
in period can both be used in determining an injection profile. Furthermore, the thermal 
data obtained when injection is resumed after a shut-in period or the data obtained from 
repeated injection and shut-in periods all can be combined to determine and/or verify an 
injection flow profile. 

[0048] An example of a modeling technique that utilizes thermal data obtained 

during a shut-in period to derive injection flow profiles is described in the following 
paragraphs. First, it should be noted that the temperature profile in the wellbore is 
affected by the fluid convection and the heat conduction between the wellbore and the 
reservoir. Because the thermal behavior of the well depends on the temperature 
distribution around the wellbore, a refined grid block scheme can be used in modeling. As 
illustrated in Figure 14, the model utilizes a grid system 132 that extends in the x, y, and z 
directions. By refining the grid size around the wellbore 30 the temperature profile can be 
stabilized. In other words, the model can utilize a grid system having a grid size selected 
such that further refinement of the individual grid sizes does not affect the temperature. 

[0049] In this type of model, it can be assumed that the temperature distribution in 

the reservoir at the shut-in has the shape of two distinctive regions, one with average 
reservoir temperature and one with the temperature of the wellbore at the shut-in. The 
temperature behavior at the wellbore can be expressed as: 




(1.27) 



The solution of this equation is illustrated in Figure 15 in terms of the dimensionless 
temperature and dimensionless time. Assuming that the temperature front is the same as 
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the saturation front, this solution can be used to relate the total injection liquid, water in 
this example, into each interval and the radius of the temperature front is as follows: 

qfa =0(1-0, (1.28) 

Or, in terms of total injection rate: 

Q inj =Af j L i t Di , (1.29) 

1=1 

This equation is used to estimate the flow profile along the wellbore with tot determined 
from Figure 15. This method relates injection rate to the temperature in an exponential 
relationship. 

[0050] Thus, a procedure for estimating the injection profile based on thermal data 

obtained during a shut-in period can be summarized as set forth in the flowchart of Figure 
16. First, the initial temperature of the reservoir, T r , is obtained, as illustrated by block 
134. The dimensionless temperature for each interval of grid system 132 is determined 
from the temperature log and the reservoir temperature: Td=(T r -Ts)/(Tr-T in j), as illustrated 
by block 136. For each unit depth of grid system 132, the dimensionless time to can be 
determined from the graph illustrated in Figure 1 5, as set forth in block 138. The 
dimensionless time for each interval along wellbore 30 is used in the equation for 
estimating the correction coefficient A, as illustrated by block 140. The flow rate into the 
reservoir for this interval is: A*toh as illustrated by block 142, and these flow rates can be 
combined to determine the overall injection flow profile. 

[0051] Accordingly, models, such as those described above, can be used to enable 

the determination of injection flow profiles in deviated wells, such as horizontal wells. 
The use of temperature sensing systems, such as distributed temperature sensors, further 
enable the desired collection of thermal data utilized by the models in deriving accurate 
injection flow profiles. 
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[0052] Although only a few embodiments of the present invention have been 

described in detail above, those of ordinary skill in the art will readily appreciate that 
many modifications are possible without materially departing from the teachings of this 
invention. Accordingly, such modifications are intended to be included within the scope 
of this invention as defined in the claims. 



21 



